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We construct the gauge invariant fermion correlator in the Schwinger model on the 

torus. At zero temperature, this correlator falls off with a rate given by the Coulomb 
£?*• energy of an infinitely heavy charge. At high temperature, the screening mass approaches 



7rT/2, and this in the presence of a mass gap. The fractional Matsubara frequency arises 
from the action of a pair of induced merons at high temperature that are localized over 
a range on the order of the meson Compton wavelength 1/m = \pKJg. We discuss the 
quenched approximation in this model, and comment on the possible relevance of some 
of these results to higher dimensions. 
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1 Introduction 

A fundamental aspect of QCD is that chiral symmetry is spontaneously broken at low 
temperature and restored at high temperature. As the low temperature phase is also 
confining, the low temperature excitations are hadronic. The nature of the excitations 
above the phase transition is still debated. The phase is believed to be screened, but 
lattice simulations indicate strong space-like correlations at all temperatures (see Ref. |1[] 
and references therein). 

Lattice simulations as well as analytical calculations have mainly focused on the be- 
havior of the hadronic correlators in the space-like (Euclidean) domain at low and high 
temperature |2|, §]. The mesonic and baryonic correlators were found to asymptote 2ttT 
and 3irT respectively, with the exception of the pion and the sigma @ . Lattice simulations 
for the quark correlator in the Landau gauge have also been performed Q. 

In non-gauge theories the fermionic propagator carries interesting information on the 
single as well as collective excitations of the system. In gauge theories the situation 
is somewhat unclear, since the naive fermion propagator is gauge dependent. In two 
dimensional QCD, 't Hooft has shown that in the axial gauge the quark propagator is 
infrared sensitive and zero in the infrared limit ||. Einhorn later noted that the infrared 
sensitivity was directly related to the gauge sensitivity of the quark propagator || . This 
feature was also noted by Casher, Kogut and Susskind in the context of the Schwinger 
model [^]. The fermion propagator has been used to probe the high temperature phase 
of QCD @,§. 

Recently, we have analyzed the mesonic correlation functions in the Schwinger model 
at finite temperature both in real and Euclidean time [[|. All correlators were found to 
fall off with the meson mass m = g/ ' ypK in the spatial direction at high temperature (for 
x 3> 1/m). Similar correlators were also investigated by Abada and Schrock [|l(| and 
Smilga |TIJ with somewhat similar conclusions. In this paper, we would like to extend the 
analysis to the gauge invariant fermion correlator in Euclidean space. By construction, 
this correlator is a combination of a fermion propagator and a path ordered exponential 
(string) — thus gauge invariant and free from the ambiguities discussed above. The gauge 
invariant correlator reduces to the ordinary fermion correlator for (time-like) axial gauges, 
barring the difficulties associated with these gauges in perturbation theory [O . Also, the 



gauge invariant correlator offers a suitable framework for probing both chiral symmetry 
breaking and confinement (screening) at zero and finite temperature through its short 
and large distance limits [0. To illustrate these points, we will use the Schwinger model 
on the torus to calculate exactly this correlator both at zero and finite temperature. At 
zero temperature, the diagonal part of the gauge invariant correlator is saturated by the 
fermion condensate at short distances, and falls off exponentially at large distances. The 
fall-off rate is related to the screening length of the attached string, a direct consequence 
of the screening character of the QED ground state. This aside, the singularities of the 
gauge invariant correlator reflect on a mass gap in the spectrum. At high temperature, 
this correlator asymptotes exp(— ttTx/2) (for x ^> 1/m) in the spatial direction and has 
a free field behavior along the temporal direction despite the fact that the spectrum 
exhibits a mass gap at all temperatures. We show that the deviation from irT in the 
spatial asymptote is due to a pair of induced merons. At high temperature, the merons 
are localized over a range on the order of the meson Compton wavelength 1/m = y/n/g- 
We discuss the quenched approximation and suggest that it cannot be applied to this 
model. The possible relevance of some of these results to higher dimensions is discussed 
in our concluding remarks. 

2 The Invariant Correlator 

Consider the Euclidean gauge invariant fermion correlator on a strip of temporal length 
j3 = 1/T, and spatial length L 

S F {x,/3) = -(q(0)e is /o^^q(x))^. (2.1) 

Although this correlator is gauge invariant, it depends on the choice of the path between 
zero and x. Below, we will always choose the shortest path on the torus staying within 
the interval [0, f3] in the time direction and [0,L] in the space direction. The expectation 
value is over the QED measure on the torus (see Sachs and Wipf [13j for more details). 
On the torus, the gauge field, A, obeys the general (Hodge) decomposition (V = (3L) [Q 

2tt . kx 1 hn h± x , 

A,(x) = -(-— + J, f ) - e„,M- (2-2) 



The first part refers to the instanton in two dimensions with topological charge k, the 
second part refers to the constant modes associated with the Polyakov line, and the final 
part is just the transverse gauge field related to the electric polarization. The calculations 
to follow are understood to be in covariant gauges [13 . 



The correlator (|2.1|) receives contributions only from the k = 0, ±1 sectors. In higher 
fluxes, the zero modes, not being absorbed by the quark fields, will make the fermion 
determinant vanish. The k = sector will contribute to the off-diagonal part (since 
the Green's function anti-commutes with 75) and the diagonal elements will be from the 
k = ±1 sectors. Inserting the explicit form of the gauge potential into ( |2.1| ) shows that 
the line integral gives a phase factor from the harmonic part (the /i's) and an integral 
from the classical background field depending on the flux factor k. The latter integral 
vanishes for paths exclusively along either the time or space direction^. The 0-dependent 
term 

igep, £ d^^-m (2.3) 

represents the quantum fluctuations around the background field and must be averaged 
in the path integral. 

The spinor structure in two dimensions allows to be factored out of the Dirac operator 

jj) = e friB<ify e !nB* (2.4) 

where Jp is the Dirac operator in the instanton background. With this in mind, the k — 1 
sector only contributes to the upper left entry of the matrix, fl2.1|) . The explicit form is 



Jd^h^\(x)MO)det\0) 1 )JD<Pe-^^ n2 - m2n ^-^ m2v e^Io XA ^ 
The zero modes are (for arbitrary positive flux k and p running from 1 to k) 

i> P ,k{x) = e- 9 ^ p , k (x) (2.6) 



-"^For other choices of paths this term would vanish upon the limit L — > 00 anyway. 



with 

U( X ) = e ^(h x°/p +hl xi/L)^ T = L jp^ 

~ 1 x° + ix 1 
h = h -p+-, and z = . 

For negative flux, —k, the sign of </> must be changed and ip(x) — > ip*(—x). The determi- 
nants of Ip have been evaluated in Sachs and Wipf |T3J]. The phase factor from the zero 
modes, U(x), cancels with the phase factor from the gauge line integral. The gaussian 
integrals can be calculated by completing squares (see Appendix A). 

After integration over <fr in both numerator and denominator, Sp,++{x, (3) reduces to 
an integration over only the harmonic part of the potential. 



_ V 2r M 2g 2 K xx -2-K/m 2 V-7T(x 1 ) 2 /V-nix°x 1 /V I 3 (x,/3) 

x I d 2 he- 27r ~ h ° T+27rfhQl $ 3 {z + hi + ir~h ) tf 3 (/n - irh ) (2.7) 



As in [13], rj refers to Dedekind's eta-function evaluated at it. The h\ integration may be 
done after expressing the theta functions as Fourier series. This gives a Kronecker delta 
that combines the two sums into one. This one sum then extends the h integration from 
the segment [0, 1] to (— oo, oo) which allows a Gaussian integration to be done. Identifying 



the fermion condensate 13 



(qq) = -^n-e ig ' K '»e-'' iK/m ' v , (2.* 



2\V\ r 2g 2 K XXr ,-2Tr/m 2 V 



the answer for the k = 1 sector is 



=i (qq)/ 



S F ,++(x,j3) = (q(x)e i9 fo A ^ q (0))';= 1 = ^^ e hW)-^? + ^?)l<>v (2 . 9 ) 

with 1% given by 



I 3 (x, 13) = 9 - £ d^d^UK^ . (2.10) 

Since OK — (□ — m 2 )~ 1 , the integration may be carried out by expanding in complete 
eigenfunctions of as done in Appendix B. The result, for L — > oo, in the temporal 



direction is (t > 0) 



h{t,P) 



irmt m 



dk 



-tVW- 



4 2 Jo [k 2 + m? 

and in the spatial direction 



7 (cosh t\Jk 2 + m 2 — 1 



/3\/fc 



2 -Urn .2 



:2.iD 



J 3 (x 1 ,/3) = -m 2 



dfc sin^^f- /3 



— coth —y/k 2 + m 2 . 

o VP + m 2 A; 2 2 



(2.12) 



Note that in the spatial direction, the high temperature limit of /3 is (x l > 0) 



7TX IT 



;(l-e" 



2/3 2m/3 

This results in a screening mass ttT/2. We want to note that this factor is not from the 
Dirac string but from the K x0 term of the integration. 

For the k = sector, the propagator is needed. In the upper right corner of the matrix, 
it is U HI 



G + -(0,x)=e 9{(f,ix) - m) U\x 



i\rj\ 3 $4(2: — H) 



3 27ri/iQ2 



(2.13) 



(3 M*)MH) 

with H = h\ — irho and G |_ is just obtained by taking the complex conjugate and the 

coordinates to their negative. Writing the determinant as 



det($ ) 



MH) 



V 



'e~ 2nTh °\ 



the h integrations may be done and the square completed in the <p integration to give 

j\ \Z e 2g 2 (K xx -K x0 ) 



S Ft+ -(x,P) 



\W l "' J 3 (x,(3)-ttz 2 /2t 



Using the form for the bosonic propagator 



g 2 K xy = 7iA m (x -y) + ln| 



V 



1 vr 



$ 1 (z 1 -z 2 )* 2V 



+ — {x'-yr + 



71 



m 2 V 



found elsewhere pL 14 , the result is 



S f 



(~ R\ - - \^//3 l#l(*)l h(x,P)-2nA m (x)-n((x°) 2 + (x 1 ) 2 )/2V-nix°x 1 /V 

\?,p)-i 2 ^^e e 



(2.14) 



In the thermodynamic limit this result can be simplified further. The phase of the $1 

function is then given by 

-bp = \Mz)\ _ | sin(7T2;) | 
6 ~ ^(z) ~ sm(nz) [ b) 

which is e(x°) for purely temporal paths and —ie{x 1 ) for purely spatial ones. For L — > oo 
A m can replaced by its Poisson resummation: 



-27rA m (x, (5) = J2 M m V( x ° ~ $ n 



x 



l\2 



The result for the lower left entry in the matrix is just the conjugate of equation ( 2.14 ). 



Putting together ( |2.9| ) and fl2.14|) we finally have for the gauge invariant fermion correlator 



in the Schwinger model at finite temperature in the infinite volume limit 



^A M*J» (( e * ° U ,-,-2^W) ( e~* 



S F (x,{3) = ^e^^ Q e ! ie )+ie-^^^ \ )) (2.16) 

where 6, the vacuum angle, has been added for completeness. For x —* ( p.!6| ) is 
consistent with the operator product expansion. The explicit form for I 3 is given in ( |2.11| ) 
and ( p,12| ) for both temporal and spatial directions. Note that we can rewrite those 



expressions in the equivalent form 
TTint f cPp 



h(t,(3) =~- J^KP 2 - - 2 )^r (*(Po) + n ) (e-^ - l) (2.17) 



Pd 
along the temporal direction, and 



h{x\(3) = +J^KP 2 - ™ 2 )^ (e(po)+np) (e mxl ~ l) (2-18) 

along the spatial direction. Here np = (e^ k +m — 1) is the Bose distribution. The 
expressions ( p . 1 7| - p7T8 ) are amenable to a spectral analysis. At zero temperature, ( 2.1 8|) 



can be rewritten in the form 

nmx 1 f (Pp 2 2\ 1Tm2 m \ ( -vnx 1 i\ fcm\ 

h = -: / -7T KP ~ m 2 )— r e(p ) (e PoX - 1 (2.19) 

as expected from 0(2) invariance. The result for zero temperature in the spatial direction 
has also been obtained in the gauge Aq = [|U|. Although the Schwinger factor is absent 



in that case it seems to us that the present approach is more transparent. 



3 Dimensional Reduction 

At zero temperature, we note that the trace of the gauge invariant correlator reduces to 
(qq)i3Cos8, which is the expected fermion condensate at finite temperature and nonzero 
vacuum angle. At large Euclidean separations, it falls off as e - m7r M/ 4 i n both the spatial 
and temporal directions. We note that 

irm g 2 f dp 1 



4 2 J 2tt p 2 + m 2 

is just the Coulomb energy of an infinitely heavy charge. The Coulomb energy follows 
from the screening of the line integral present in the gauge invariant correlator in the 
QED vacuum. It is a purely "kinematical" term. This aside, the expressions (|2.17|) and 
( P-19Q show that the singularities of the gauge invariant correlator are related to the 
mass gap m with a form factor Trm 2 /p 2 . The off-diagonal part of the gauge invariant 
correlator reduces to the free fermion propagator at short distances i'y ■ x/\x\ 2 . The large 
distance behavior is similar to the diagonal part. Our result ( |2.16| ) at zero temperature 
is in disagreement with the brief analysis by Sachs and Wipf of a related correlator (|13|], 
section 6). 

At finite temperature, we note that the diagonal part of the gauge invariant cor- 
relator reduces to the temperature dependent condensate at zero separation, while the 
off-diagonal part gives the usual free fermion propagator. At high temperature, the gauge 
invariant correlator along the temporal direction exhibits a free field behavior. Along 
the spatial directions, the gauge invariant correlator falls off at a rate which is given by 
e -TrT\x\/2 d eS pjt e the f ac t that the spectrum has a mass gap m = gj \[k- 

To understand this behavior, we note that at high temperature the model dimension- 
ally reduces to one dimension since (u n = 2irnT). 

<f>(t, x)=J2 e~ iWnt 0n(x) ~ o (x) (3.2) 

n 

Therefore, A\ ~ and 

A (x)~ ^h-^Mx)- (3-3) 

With this in mind, the diagonal part of the gauge invariant correlator reduces to 

8 



(q(0,x)e ia f: dx ' A ^> x \(0,0)) p ~ (q(0, x) q(0, ())>,. (3.4) 

Because of the trace over Dirac indices, only the zero modes contribute to (|3.4j) . In the 
dimensionally reduced theory, the expectation value becomes 

(q(0, x) q(0, 0))^ ~ ( e -s«<o(-)+<MO))^ ( 35 ) 

The expectation value is with respect to the dimensionally reduced action 

S 0+1 = pjdx^ (V 4 - m 2 V 2 ) 0. (3.6) 

The integral being quadratic in can be performed at once using the saddle point equa- 
tion. 

(V 4 - m 2 V 2 ) <p 0t t(z) = -gT (S(z - x) + 8{z - 0)) (3.7) 

The solution is an induced instanton (first introduced in the context of mesonic correlation 
functions [|HJ). 

TFT / 1 . , , , 1 



<f>oJz) = — — e" m|z| + \z\ + — e-™!*-*! + U - xl (3.8) 

2g \m m J 

Inserting ( |3.8|) into (|3.5|), (|3.6|), and using (|3.7|) yields (|x| — >■ oo) 

(q(0, x) q(0, 0))^ ~ e -^(W+£+£ e - ra| ") ~ e^'^/" 1 . (3.9) 

For the off-diagonal part the induced instanton gives rise to a factor exp(7rTx/2) which to- 
gether with a factor exp(— nTx) from the free propagator results in the same asymptotics. 
As final answer we find 



(q(0, 0) q(0, x)) p ~ _ e -*m/2-,T /m f ii _ !2_ ) (3.10) 



2 2(3, 
at high temperature. 

We remark that for x ^> 1/m the induced instanton ( [3.8|) is a linear superposition of 

two merons. Indeed, the classical field 



gives rise to the following potential (ho = —l/4g) 

M*) = ^ s s n W ( e_mN - - 2J (3 ' 12) 

with A (— oo) = and A (+oo) = —iiT/g. The topological charge carried by ( |3.11| ) is 

* = g £f dxE ( x ) = -4a {M+oo) - A ° ( " oo)) = \ (3 - 13) 

one-half of the instanton charge. In the Schwinger model the merons are induced (refer 



to (3/7 )) and localized over the meson Compton wavelength 1/m since the electric field 



E(x) ~ Te m \ x \. When x becomes on the order of 1/m, the two merons merge into a single 



instanton. The solution for x = was studied in |ll[; its classical action is responsible 
for the temperature dependence of the condensate. 



Mesonic correlation functions can also be analyzed along these lines JO]. In that 



case the action of the classical solution cancels the exp(— 2irTx) asymptotics of the free 
propagator resulting in an exponential fall off determined by the meson mass. In the 
present case both the classical solution and its action are a factor two smaller (as compared 
to the mesonic action in the k = 2 sector) and in the k = ±1 sectors there is no additional 
exp(— tcTx) dependence resulting in the asymptotic screening mass nT/2. In the k = 
sector the action of the two meron configuration gives rise to a factor exp(7rTa;/2) but in 
this case a factor exp(— 7rTx) comes from the free propagator. Merons at high temperature 
therefore give rise to fractional Matsubara frequencies in the Schwinger model. 

4 The Quenched Approximation 

The quenched approximation for the Schwinger model has been discussed extensively in 



the literature |17| and a number of contradictory statements have been made. Part of the 
confusion will be discussed below. Before this, however, we would like to stress the fact 
that in the quenched approximation and at zero temperature there are no gauge degrees 
of freedom. The fermionic ground state is trivial, and as such (qq) should be zero. In the 
absence of the Dirac sea, there is no (axial) anomaly in the quenched approximation in any 
dimension. At finite temperature, and in the Euclidean formulation, the Polyakov lines 
play the role of non-local gauge degrees of freedom. Dimensional reduction arguments 



10 



suggest that the high temperature phase in the quenched approximation is still trivial. 

Having said this, we may formally calculate (qq) for fixed fermionic mass ii, number 
of flavors Np and two volume V, and consider the limits /i — > 0, Nf — > and V — > oo. 
Some of these calculations have already appeared in the literature, with some confusion 
related to the order of the limits. In general, each pair of these limits does not commute. 
Let us first discuss the mass (//) dependence. The important parameter is £ = iiV (qq). 
For (<1, the mass is much smaller than the smallest eigenvalue and, as a consequence, 
only the sectors with topological charge equal to ±1 contribute to the condensate. For 
£ ^> 1 the mass is much larger than the smallest eigenvalue, and the condensate receives 
contributions from many different topological sectors |IS) . Remarkably, in both cases the 



value of the condensate is the same. To be more precise, if £ <C 1, the primed sum in the 
expression for the condensate 

is of order jjl and vanishes with respect to the contribution from the zero eigenvalue. For 
£ ^> 1 the second term becomes dominant. 

In the quenched approximation, the gauge invariant correlator receives contributions 
from all /c-sectors since the restriction from the zero modes in the integral is no longer 
present. The same effect occurs by taking a sufficiently large quark mass, and we expect 
that in this limit the quenched approximation is valid. In the opposite case, £ <C 1, the 
condensate diverges in the chiral limit. This is clear from (|4. 1|) because the factor iii 
from the determinant is absent in this case. The extra weight in the bosonic integration 
|m 2 0D0, which originally was from the fermionic determinant [I3| is no longer there. 



This results in a different bosonic Green's function. With this in mind, the fermionic 
propagator reads _ 

q(0)qW = E *^M (4.2) 

with fx being a regulator for the zero modes which will be taken to zero after L — > oo. In 
the Schwinger model and also for QCD with one flavor the result for the condensate is 
independent of the order of the limits \x — > and V — ► oo. Assuming that this is also the 
case for the Schwinger model in the quenched approximation, we choose \x <^i 1/yV so 
that only the k — ±1 sectors contribute to the condensate in the quenched approximation. 

11 



Therefore, for (<1, the fermion condensate becomes 

2i Ek>oke- 27Tk2 / m2v 



{qq)Q ~ W j: k>0 e-^/^y ■ (4 - 3) 

In the limit V — ► oo this sum can be calculated: 

. iy/2 m ,. ,, 

(qq)g = iff (4-4) 

in agreement with an analysis using the proper time representation of the fermion propa- 
gator [19|. The fact that the condensate diverges for \x — ► seems to be in contradiction 



with arguments presented at the beginning of this section. However, the condensate as 
calculated in (4.3) follows from the trace over all fermionic states and cannot be identi- 
fied with the condensate that would be obtained in the absence of the Dirac sea (which 
is zero). 

The gauge invariant fermion correlator cannot be calculated exactly in the quenched 
approximation. The Wilson line - - being different in the quenched approximation - 
shows an area law (quenched) as opposed to a perimeter law (unquenched) . As such, 
it would be interesting to evaluate the ratio of the invariant correlator to the fermion 
condensate in the quenched approximation. 

5 Conclusion 

We have explicitly constructed the gauge invariant fermion correlator for the Schwinger 
model on the torus. We have shown that at zero temperature, the result is 0(2) invariant 
and asymptotes e~ crp ' :E ' at large separations with o> = ixm/A being the screening mass 
associated with the line integral. This screening mass is purely kinematical, and relates 
directly to the perimeter law of the Wilson loop at zero temperature 

where P is the perimeter spanned by the large loop C. At short distances, the diagonal 
part of the gauge invariant correlator reduces to the fermion condensate, while the off- 
diagonal part reduces to the free fermion propagator. After proper subtractions, the 
gauge invariant correlator displays singularities that are related to the mass gap in the 
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theory. These singularities may be related to heavy-light systems, with the string playing 
the role of an infinitely heavy charge. At finite temperature, we have shown that the 
diagonal part of the gauge invariant correlator at short distances reduces to the fermion 
condensate and asymptotes exp(— irTx/2) at large spatial separations, even though the 
spectrum exhibits a mass gap at all temperatures. The fractional Matsubara frequency 
results from the action of localized merons at high temperature. We have shown that the 
quenched approximation does not work in the Schwinger model. These results suggest 
that this is also the case for QCD with massless quarks. However, it may very well be that 
for nonzero quark masses the effect of the fermion determinant can be ignored for most 
observables. The relation of these results to the gauge invariant correlator in real-time 



will be discussed elsewhere |21J . We briefly note that the gauge invariant correlator is not 
a periodic function in the temporal direction. 

Are the above results relevant for four dimensional QCD? We do not know. We suspect 
however, that the gauge invariant correlator for QCD with one flavor should reflect on 
the fermion condensate along the diagonal. It should asymptotically approach e~ <Jp ' a; ' at 
zero temperature in all directions with <jp as the perimeter-law coefficient in the Wilson 
loop. It would be interesting to check this point by lattice QCD simulations and relate the 
results to the (renormalized) Coulomb energy in four dimensions. This kinematical term 
aside, we suspect that the rest of the correlator should fall off at a rate determined by 
the heavy-light bound states of QCD. At finite temperature, the diagonal part should be 
proportional to the finite temperature fermion condensate. In particular, it should vanish 
in the chirally symmetric phase. Lattice simulations of the mesonic screening masses 
together with dimensional reduction arguments [I], |3| suggest that at high temperature 
the screening masses asymptote 2iiT. Repeating the dimensional reduction arguments, 
and barring induced merons at high temperature, we would conclude that the invariant 
fermion correlator asymptotes irT and not a fraction of iiT in QCD. We recall, however, 
that merons in QCD have been proposed as candidates for understanding confinement 
at zero temperature ||21|| . Their introduction at zero temperature is, however, ad hoc, as 



they carry infinite action before smearing. It would be interesting to see whether they 
could be induced and localized (have a finite action) in QCD at high temperature. Lattice 
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simulations in these directions could be helpful. 
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Appendix A 

We detail in this appendix the calculations that led to (2/7) in the text. The factors 
associated with the exact Dirac operator in ( |2.4|) , along with the line integral ( |2.3|) , act 
as effective sources in the action for the <fi field. 

5 eff [0] = y>y{^(D 2 -m 2 D)0 + 0J} 

= - f d 2 y{4>(n 2 - m 2 n)0 - J(n 2 -m 2 n)- l j}, (A.l) 



with 



1 f x d 

-J{y) = ej 2 (y - x) + e £ 2 (?/) + ie^ / d^——5 2 {y - f ). 

o jo at 



-Jo ^.9^ 

Here, the e^'s are either plus or minus depending on in which fc-sector the correlator is 

evaluated. They have the same sign in the k = ±1 sectors and opposite signs for the 

propagator in the k = sector ( |2.13| ). The fermions contribute strongly at the endpoints 

(as delta functions) but the gauge field's contribution is smeared out along the transverse 

direction of the path. 

The source-source term in equation (|A.1|) given in terms of the bosonic Green's func- 



tion, K = (D 2 — m 2 D) 1 , and including the extra minus sign in the exponent is 

g 2 (K xx + e x e K x0 ) + ie x l x x + %e Q l\ + I 2 (A.2) 



with 



Jo 



10 

,2 



9 f x f 

h = —z^tpa / / d^dr] p d^d^ K iv . 

Li J •/ U 

The first term in equation ( |A.2[) is just the self-energy contribution from the fermions 



along the path. The vanishing of I\ follows from the mode expansion of the Green's 
function. Due to the identity e pu e pa = 5 w 5 va — S^ a 5 up , the self-energy of the gauge field, 
I2, may be rewritten as follows 

9 

h = g\K xx - K x0 ) + 9 - f d^d^UK^ = g 2 (K xx - K x0 ) + J 3 (x, /?). 

z Jo 

giving finally in our case for the bosonic contribution in the exponent 

15 



2g 2 (K xx + K x0 ) + I 3 (x,(3) 

in the k — ±1 sectors and 

2g 2 (K xx -K x0 ) + I 3 (x,P) 

in the k = sector. 

Appendix B 

In this appendix, we detail the calculation that leads to the explicit form of the I 3 {x, /3) 
integral both along the temporal and spatial directions. The double integral of the bosonic 
Green's function needed in Appendix A is 

7 2 rt 



h{t,f3) = 9 - j dt'dt"nK(t'-t",x)\ x= 

2 Jo 



in the temporal direction. Since K = (D 2 — m 2 D) _1 , an additional □ in the numerator 
leaves only the massive part of the propagator. Using a complete set of states of (f> this 
may be rewritten as 

nm 2 ft . t e 2mno(t'-t")/l3 



h(t,(3) = -^- dt'dt" £ 



2V J ° J ° (no.ntjJco.O) (2™o//?) + (2^/ L)' + m 2 

7rm 2 t 2 *-^ 1 



E 



2V ^ (27rn/Lr + m 2 



ra^O 



7r m(3 , 2 ^ >r-^ sin 2 Tint/ (3 

i_ v o ' 2.^1 2.^1 



It 2 > ni ^ oon ^ (7rno) 2 ((7rn ) 2 + (7rn 1 /r) 2 + (m/5/2) 2 ) 

The no sum in the second term may be done by integrating both sides of 

cos27m(/> 1 cosha(l — 2|0|) 1 



y , 

~^ (7m) 2 + a 2 a sinha a 2 

with respect to twice. The result in the temporal direction for finite L is 

Tit 2 nm\t\ , mi 7rm 2 J5, coth f£(n) (cosht£(n) — 1) — sinh |£|£(n) , N 

w-t^t-tS ?w (B1 ' 



with £(n) = J(2irn/L) 2 + m 2 . 
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Taking the limit L — > oo, the sum turns into an integral and gives 

r , „ s Trm\t\ m 2 r°° dk r . 3 r-r, ~, . rn\ n ,-, 

I 3 (t,P) = -^ / r coth -Vfc 2 + m 2 (cosht\/fc 2 + m 2 - 1) 



nm\t\ 


m 2 f°° dk 
2 Jo (fc 2 +m 2 )3 


coth— -\ 


4 




— sinh t a/A; 2 + m 2 



The part in the brackets may be rewritten as 

e -\t\Vk*+m* _ l + — _^_ ( cos h tVF + m 2 - 1) 

showing the Bose-Einstein number occupation factor for finite temperature with a form 
factor. 

For the spatial result, all that is needed is to exchange x with t and L with /3 in 



equation B.l. Taking L — * oo in that result produces 



/ 3 (^,/3) = ^coth-^ + — - J2 



4 2 2/3 n fr' 00 ((2vmT) 2 + m 2 ) 2 

2 r 00 rfA: sin 2 ^h /5 / — - 

io Vfc 2 + m 2 k 2 2 

In the first expression T — 1/3 was used to show the Matsubara form. The second 
expression may be obtained by expanding the hyperbolic cotangent in the first expression 
as a series, combining the two series, and noting that 

r 00 dk sin 2 '' "'" 



r°° dk sin -y -na n 
Jo k 2 (k 2 + b 2 ) = ^ _ 4& 2 ^ ~ 6 



/o k 2 (k 2 + b 2 ) b 2 Ab 2 
to produce the integral. 
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